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ABSTRACT 

Infrared Dark Clouds are expected to harbor sources in different, very young 
evolutionary stages. To better characterize these differences, we observed a sam- 
ple of 43 massive Infrared Dark Clouds, originally selected as candidate high-mass 
starless cores, with the IRAM 30 m telescope covering spectral line tracers of low- 
density gas, high-density gas, molecular outflows/jets and temperature effects. 
The SiO(2-l) observations reveal detections toward 18 sources. Assuming that 
SiO is exclusively produced by sputtering from dust grains, this implies that at 
least in 40% of this sample star formation is on-going. A broad range of SiO line- 
widths is observed (between 2.2 and 65kms _1 ), and we discuss potential origins 
for this velocity spread. While the low-density tracers 12 CO(2-l) and 13 CO(1-0) 
are detected in several velocity components, the high-density tracer H 13 CO + (l- 
0) generally shows only a single velocity component and is hence well suited for 
kinematic distance estimates of IRDCs. Furthermore, the H 13 CO + line-width is 
on average 1.5 times larger than that of previously observed NHs(l,l). This is 
indicative of more motion at the denser core centers, either due to turbulence or 
beginning star formation activity. In addition, we detect CH 3 CN toward only six 
sources whereas CH3OH is observed toward approximately 40% of the sample. 
Estimates of the CH3CN and CH3OH abundances are low with average values 
of 1.2 x 10~ 10 and 4.3 x 10 -10 , respectively. These results are consistent with 
chemical models at the earliest evolutionary stages of high-mass star formation. 
Furthermore, the CH3OH abundances compare well to recently reported values 
for low-mass starless cores. 

Subject headings: stars: formation - stars: early-type - ISM: jets and outflows 
- ISM: abundances - ISM: clouds - ISM: kinematics and dynamics 
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Introduction 



While observational research in high-mas star formation over the last few decades has 
focused on relatively evolved regions like Ultracompact Hll regions ( UCH l ls) and High-Mass 
Protostellar Object s (HMPOs) (e.g., see reviews by Kurtz et al.l I200CH : IChurchwelll 12002 



Beuther et al.l 120071 ). the earliest stages of massive star formation have become regularly 
accessible only through the advent of large-scale mid-infrared surveys of the Galactic plane 
by satellite missions like the Infrared Space Observatory (ISO), the Midcourse Space Ex- 
periment (MSX) and most recently Spitzer. The basic observational characteristics of the 
earliest stages of (massive) star formation, prior to the formation of any embedded heat- 
ing source, should be that they are strong cold dust emitters at (sub)mm wavelengths, and 
weak or non-detections in the mid-infrared because they have not yet heated a warm dust 
cocoon. In the low- mass regime, so-called Infrared Dark Clouds (IRDCs) have been studie d 



regularly since a few years (e.g., lEgan et al.l 11998c ICarey et al.ll2000c iBacmann et al.l 120001 ) 



Sridharan et al. 


2005 


Hill el 


: al. 


2005; 


Klein et al. 


2006; 


Rathborne et al. 


2006; 


Pillai et al. 


2006h. 



2005; Simon et al. 2006: Beltran et al. 



(e-g-,, 



Massive IRDCs (in short IRDCs from now) are not a well defined class, but they likely 
harbor various evolutionary stages, from genuine High-Mass Starless Cores (HMSCs) via 
High-Mas s Cores harboring a ccreting Low/Intermediate-Mass Protostars to the youngest 
HMPOs (IBeuther et al.l 120071 ). While the first stage provides good targets to study the 
initial conditions of massive star formation prior to cloud collapse, the other stages are 
important to understand the early evolution of massive star-forming clum ps. Here, we are 



target ing a sample of massive IRDCs first identified as candidate HMSCs by lSridharan et al. 



( 120051 ) to differentiate between the different evolutionary stages existing within this sample. 
To facilitate this, we observed the target sample with the IRAM 30 m telescope in different 
spectral lines sensitive to various physical processes like molecular outflows (e.g., SiO, CO), 
dense gas (e.g., H 13 CO+, CH 3 CN, CH 3 OH) or hot core conditions (CH 3 CN, CH 3 OH). Table 
[U lists the observed spectral lines, the upper level excitation energies -E U p P er and the critical 
densities n cri t- 



2. Observations and Data Reduction 



We observed a sample of 43 IRDCs selected from lSridharan et al.l (120051 1 in August 2005 
at the IRAM 30 m telescope on Pico Veleta with four receivers simultaneously. The flexible 
multi-receiver setup of the 30 m telescope allowed us to observe usually a few lines within 
the same spectral setup of each receiver ((1) H 13 CO+(1-0) & SiO(2-l), (2) 13 CO(1-0) and 
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CH 3 CN, (3) 12 CO(2-l), (4) CH 3 OH). TableHlists the corresponding frequencies, the spectral 
resolution and the system temperatures for each line. The observations were conducted in a 
position-switching mode with OF F positions carefully chosen from the Massachusetts-Stony 
Brook Galactic plane CO survey fjSanders et al.lll986l ). usually approximately of the order 
1000" offset from the source positions. 



Table 1: Observing parameters 
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a Average system temperature and lcrrms values. 

b Calculated from A = 0.3A7p 3 nM 2 (Amn in units of 100 yum and fi 
7 (CH 3 CN(6 - 5)) 20K from lPei fe Zem3 (|l995bl Jah. 
c ^crit = A/'j with the Einstein coeffi cient A and the col lisional rate 7. Except for CH 3 CN, 
A and 7 were taken from LAMBDA JSchoier et al.ll2005h at T = 20 K. 
d No data for this transition in LAMBDA. 



3. Results and discussion 

3.1. 12 CO(2-l), 13 CO(1-0) and H 13 CO+(1-0) 

Figures [H [2] and E] show the 12 CO(2-l), 13 CO(1-0) and H 13 CO + (1-0) spectra for the 
whole source sample. Except of IRDC 18431-0312-4, where the H 13 CO + (1-0) line remained 
undetected, toward all other sources the three lines were detected. The most obvious dif- 
ference between 12 CO/ 13 CO on the one side and H 13 CO + on the other side is that the 
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CO isotopologues are detected at many different velocities whereas H 13 CO + usually only 
shows a single spectral peak (except of a few double-peaked structures listed in Table [2]). 
This difference can be explained straightforwardly by the location of the sources within the 
Galactic plane because CO traces more or less all molecular clouds along the line of sight, 
hence contributions from various Galactic spiral arms, whereas H 13 CO + has significantly 
higher critical densities (Table [1]) and thus traces only the dense cores we are interested in. 
Therefore, if one wants to know the velocity of a specific region CO observations are mostly 
not useful and one has to refer to high density tracers. The peak velocities v of H 13 CO + 
(Table [2]) corr espond well to the pre viously derived velocities from NH3 observations of the 
same sample ( Sridharan et alJbood ). Identifying the correct associated velocity component 
to derive accurate velocities i s an important s tep to get distance estimates from the Galactic 
rotation curve as reported in lBrand fc Blitzl (119931 ) and ISridharan et al.l (120051 ). 



Table [2] also lists the Full Width Half Maximum (FWHM) line-widths Av (H 13 CO + ) of 
the observed H 13 CO + lines from Gaussian fits, the average H 13 CO + line-width is 2.6 kms -1 . 
FigureHJshows a comparison of the H 1 3 CO + line-width with the previously derived line- width 
from the NH 3 (1,1) observations from (ISridharan et al.ll2005l ). We find that the H 13 CO + line- 
width nearly always exceeds the line-width from the NH 3 (1,1) with an average line- width 
ratio Aij ^(nh^(ii)) ^ = This is interesting since the critical densities f or both molecules 



differ by about 2 orders of magnitude (ra crit (NH 3 (l, 1)) ~ 2 x 10 3 cm 3 , e.g. IZhou et al.lll989 
n crit (H 13 CO + (l - 0)) ~ 1.7 x 10 5 cm" 3 , Tabled]), whereas the thermal line-widths (Av th = 
81n2-^- with the molecular mass m mo \) are similar (AuthfNlWl, l))i5K ~ 0.2 kms~ 

^rriol 



and Av th (H 13 CO+(l - 0))i 5K ~ O.lSkms" 1 ), Aw th (H 13 CO+(l - 0)) being even smaller than 
Af t h(NH 3 (l, 1)) because of its larger molecular weight. Therefore, differences in the non- 
thermal motions have to be responsible for the line-width variations between both species. 
Similar results we re obtained for low-mass cores in a comparison of CS with NH 3 line- widths 
(IZhou et al.lll989h . 



There exist several potential ways to explain such a line- width differenc e. For example, 
there is the well-known line- width-size relation Av oc R 5 ( iLarsonl Il98ll ). This relation 
would imply that our target regions would be on average 2.3 times larger in H 13 CO + than 
they are in NH 3 . Although we have no maps of the re gions , this appears very unli kely. A 
different possibility was discussed by IZhou et al. Jl989h and lMvers fc Benson! Jl983h : If the 
stability of a core is determined by the entire line-broadening, for a gravitationally stable 
core the critical density should relate to the line-width and size of a molecular core as 



n crit oc (Av/Rf 



Av oc R^n crit . 



To fulfill this criterion the H 13 CO + source size would need to be ~7 times smaller than 
that of NH 3 which again appears to be very unlikely. In this scenario, one could then 
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conclude that molecular clumps are not gravitationally stable but maybe collapsing already. 
However, overall collapse is not sure at all from such an anal ysis, and the clum ps could also 



be fragmenting or only partially collapsing (see discussion in lZhou et al.lll989l ). 



As discussed below in §3.21 a significant fraction of sources exhibits broad line wings in 
the jet/outflow tracer SiO, which is indicative of early outflow and hence early star formation 
activity. Since H 13 CO + has a larger critical density it traces denser regions of the gas clumps. 
With ongoing star formation activity, one expects additional line broadening effects from 
various physical processes like outflows, infall and potential rotation of embedded accretion 
disks, all spatially unresolved by the current observations. Most of these effects are likely 
contributing more strongly in the dense gas components and hence could explain the excess 
in H 13 CO+ line-width. 



3.2. SiO(2-l) 

One of the aims of this study was to identify how many of the sources would show clear 
SiO emission caused by molecular outflow activity. SiO is believed to be produced in the gas- 
phase after sputtering Si-bearing speci es from the dust gr ain surfaces when the outflow/jet 



impinges on the ambient gas and dust (jSchilke et al.lll997l ). SiO is usually barely detectable 



in quiescent clouds. Figure H] presents the 18 SiO(2-l) detections toward our source sample 
of 43 regions. If this source sample is representative for typical Infrared Dark Clouds, this 
implies that in at least 40% of the IRDCs star formation has already started (absence of 
SiO emission does not necessarily imply absence of star formation activity). However, it is 
not obvious whether this sample is representative for all IRDCs because it was selected to 
be close to sites of already ongoing massive star formation. 

Table [2] also lists the width-to-zero-velocity Avq of the observed SiO(2-l) lines. If one 
divides the observed SiO spectra by their width-to-zero-velocity Av into 3 regimes (a) low 
velocity At> < lOkms -1 , (b) intermediate velocity lOkms -1 < Av Q < 20kms _1 and (c) high 
velocity Av o > 20 km s _1 , one finds six sources in the low- velocity regime, five at intermediate 
velocities and seven in the high-velocity regime. While the sources at intermediate to high 
velocities obviously show gas at relatively high velocities with respect to the ambient gas, 
this is less clear for the sources in the low- velocity regime: those sources with Avq > 5 km s" 1 
probably also have outflow contributions, but the SiO (2-1) spectra toward two sources with 
Av < 2.5 km s" 1 (IRDC 18223-6 & IRDC 18530-2) hardly look like typical ou tflow spectra. 



A similar spread in SiO v elocities has also previously been observed by, .e.g., iLefloch et al. 



19981 ; ICodella et al.lll999l . What are possible reasons for the different observed SiO Avq c : 
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The observed velocity spread could partly be produced by the angles between the out- 
flows and the plane of the sky since we only observe line-of-sight velocities. However, as- 
suming a random distribution of outflow inclination angles 9 (9 = is an outflow along the 
line of sight) and similar unprojected outflow velocities, one would expect fewer outflows at 



high velocities than at low velocities (J 45 sin(9)d9 ~ 0.29 vs. sin(9)d9 ~ 0.71). This 
is in contrast to our relatively even distribution of Avq from low to high values. Hence, 
inclination effects are unlikley to explain the observations. 

Outflows/jets from massive p rotostars have intri nsically higher velocities than those of 
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their low-mass counterparts (e.g., iRicher et al.l 120001). In additio n, more luminous sources 
usually exhibit stronger SiO emission (e.g., ICodella et al.l Il999l ). Therefore, the sources 



with smaller line-width could harbor less massive protostars that may or may not form a 
massive star at the end of their evolution. Furthermore, all objects with low SiO velocity 
spread are also the weakest or non-detections in CH3CN and CH3OH (see §3.3p which is 
also expected in the framework of embeedded protostars with different masses. Thus, differ- 
ent masses/luminosities of the embedded outflow-driving sources can potentially cause the 
observed SiO line-width difference. 

Since IRDCs are believed to represent the earliest evolutionary stages of massive (and 
low-mass) star formation, we expect that in some of them low- to intermediate-mass proto- 
stars may be embedded that are likely destined t o grow, accrete and become high-mass stars 
at th e end of their evolution (e.g., IRDC 18223-3. iBeuther et al.ll2005bl ; iBeuther fc Steinacker 
20071 ). However, a priori we cannot differentiate whether an object will become massive at 
the end of its evolution or whether it will remain a low-mass protostar. In this scenario, one 
can also speculate about an evolutionary sequence where the growing protostars accelerates 
the gas to higher velocities in the course of its early evolution, i.e., the smaller line-widths 
may correspond to the still less massive and younger objects whereas the broad line-width 
sources may have already accreted more mass and could hence be a little bit more evolved. 
This picture is consistent with the above assertion that the line-width distribution can de- 
pend on the mass and luminosity of the central object, just with the additional assumption 
that the lower-mass objects may still be able to become massive stars during their ongoing 
evolution. This scenario may appear counter-intuitive compared to typical low-mass outflows 
where the strongest emission is observed from class sources and whe re the SiO emission gets 
weaker formore evolved class I and II sources (e.g.. lGibb et al.ll2004j ). However, even class 
sources have typically already ages of th e order 10 4 yrs, and IRDCs l i ke those presented her e 
are likely still significantly younger (e.g.. IBeuther fc Steinackerl 120071 ; iKrumholz et al.ll2007l ). 
Furthermore, this low-mass classification scheme is unlikely to be well applicable to high- 
mass star formation because, among other reason s, massive stars reach the main sequenc e 
while accretion processes can still be ongoing (e.g.. lSchaller et al.lll992l ; iMcKee &: Tanll2003l ). 
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Therefore, our picture applies only to the earliest evolutionary stages of strongly accreting 
protostars destined to become massive stars at the end of their evolution. In principle, the 
non-detection of SiO line- wings could also be caused by poorer signal-to-noise ratios or larger 
distances to the respective region. However, both scenarios are unlikely the case here because 
we have approximately the same rms for all sources, and the line-width differences are found 
also toward sources at the same distances (e.g., IRDC 18223-1243-3 & IRDC 18223-1243-4 
versus IRDC 18223-1243-6). 

Do the Avq < 2.5 kms -1 fit into any of these pictures? Since more than 50% of the target 
sources exhibit no SiO emission at all down to our s ensitivity limi t, it is intere sting that two 



object s show these extremely narrow line-widths. iLefloch et al.l (119981 ) and ICodella et al. 



(119991 ) try to explain similar observations by outflow/dense clump interactions and in an 
outflow evolutionary model. In the latter scenario, the highest velocity components are 
found at the earliest evolutionary class stages and the smaller line-width components 
could be remnants from earlier interactions of molecular outflows with the ambient gas 
(jCodella et al.l Il999l ). They estimated the corresponding time-scale to decrease the SiO 
emission and the velocities to be of the order 10 4 yrs. We consider this scenario as unlikely 
for our sample b ecause massive star f ormation proceeds on relatively short time-scales (of the 
order 10 5 years, iMcKee fc Tanll2002l ). and after 10 4 yrs the regions are expected to be strong 
en ough to be detetable a t near- and mid-infrar ed wavelengths (e.g., the samples presented 
by lMolinari et al.lll996l or ISridharan et alihood ) which is in contrast to our sample selection 
as Infrared Dark Clouds. Another possibility is to speculate whether these objects could 
represent the earliest observed stage of star formation within our source sample: In this 
picture, the Si-bearing species would just be released from the grains by shock interaction 
and SiO would only have begun to form in the gas phase, but the young jet/outflow would 
not have had the time yet to accelerate large amounts of gas to detectable higher velocities. 
To confirm or discard this picture additional indicators of star format ion activity would 



need to be observed, e.g., infall of gas via asymm etric line profiles (e.g., iMyers et al.lll996l ) 



or shocked H2 emission either in t he K-band (e.g.. lDavis et al.ll2004j ) or in the IRAC 4.5 /im 



band (e.g.. iRathborne et al.ll2005l ). 



3.3. CH 3 CN(6 fe - 5 fc ) and CH 3 OH(5-4) 

In addition to the previously discussed molecules, the spectral setup allowed to observe 
also the typical hot-core tracer CH3CN as well as a series of four CH3OH lines. Although 
the critical densities and excitation parameters of both molecules do not differ that much 
(see Table [TJ, CH 3 OH typically has one to two orders of magnitude larger gas column 
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densi ties and is hence easier to observe (e.g.. Ivan Dishoeck Blakdll998l ; iNomura & Millar 
20041 ). Furthermore, CH 3 OH is known to be present in massive star formation right from the 
beginning, hence it is a so-called parent molecule, whereas CH 3 CN is a daughter molecule 
which needs t ime to form through c hemical processes during the massive star formation 



process (e.g., INomura fc Millarl 120041 ) . Therefore toward a sample of young IRDCs, one 
would expect CH 3 CN to be either undetectable or only weakly detected in a few sources, 
whereas CH 3 OH should be detectable toward a larger fraction of the sample. 

These expectations are largely confirmed by our observations (Figs. [6] & CJ), we de- 
tect CH 3 CN only toward 6 sources (~14%) whereas CH 3 OH is observed toward 17 regions 
(~40%). Furthermore, if detected, CH 3 CN(6fc — 5&) is observed only in its two lowest k- 
components with excitation temperatures < 30 K, the k = 2 transition is detected in no 
source at all. Similarly, for CH 3 OH, we detected toward the 17 regions the two lower ex- 
citation lines, whereas the two higher excitation lines are detected only toward a small 
sub-sample of 5 regions. This is additional confirmation of the youth and low average tem- 
peratures of the selected source sample. Table [3] lists the line-widths Auk=o and Az/fc =i of 
the CH 3 CN k — 0,1 components from Gaussian fits to each line. While we cannot fit one 
source because of multiple components and in one other source Av k=0 and Au k=1 are of 
approximately equal width, for the four other sources we find consistently larger line-widths 
in the k = 1 component compared to the k = component. This is consistent with already 
ongoing star formation activity in these regions, where the warmer more central gas is sub- 
ject to more non-thermal motion, either of turbulent nature or due to infall, rotation and/or 
outflow processes. 

One can try to use the observ ed lines also to ge t more quantitative temperature and col- 
umn density estimates. Following lLoren fc Mundyl (119841 ). assuming optically thin emission 
and that a single temperature characterizes the populations within a fc-ladder and between 
different /c-ladders, the ratios of the intensities of different k lines can be used to estimate 
this temperature. Specifically, for the J = 6 — 5 line and k = and k = 1 components, we 
have 



T = 7.2/ln(I /h) 

with Ik the corresponding observed peak intensities. For IRDC 18102— 1800-1 and 
IRDC 20081+2720-2, the relative line intensities are consistent with optically thin emis- 
sion and we estimate the temperatures to be 16 ± 16 and 16 ± 13 K, respectively. For the 
other sources, while they are unlikely to be optically thick, the intensities are inconsistent 
with optically thin emission, so we do not estimate temperatures. It is interesting that the 
estimates, where possible, although not very accurate, turn out to be low and similar to 
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estimates from NH3 observations (ISridharan et al.ll2005l ). We cannot use the two mainly 
detected CH 3 OH lines for temperature estimates because they are of different symmetry (A- 
and E-CH3OH) and can be considered for radiative transport purposes as distinct molecules. 

However, we can calculate the column densities in the detected CH3OH and CH3CN 



lines, and then derive the total molecular column densities following iRohlfs fc Wi 



As pa rtition functions we used the values listed in the JPL catalog at T = 18.75 K (jPoynter fc Pickett 



son J 2006). 



19851 ). Furthermore, we are interested in CH 3 OH and CH 3 CN molecular abun dances. To es- 



timate these, we ex tracted the 1 .2 mm pea k fluxes from th e image s presented in lBeuther et al. 



(b002j fl Following Eildebrandl Jl983h and lBeuther et al.l fl2005ah . we then calculated the H 2 
column densities at T = 19 K (corresponding to the partition functions for the molecular 
emission) assuming optically thin emission and a dust opacity index (3 = 2. The spatial 
resolution of the 1.2 mm and CH 3 OH data is approximately the same, for the comparison 
with the lower-resolution CH3CN data, we smoothed the 1.2 mm images also to the corre- 
sponding spatial resolution of 23". The 1.2 mm peak fluxes and H2 column densities (both 
at 1 1" resolution) are listed in Table [2j The different H2 column densities were then used to 
estimate the CH3OH and CH3CN abundances. The calculated molecular column densities 
are estimated to be correct within a factor 5, the accuracy of the estimated abundances 
is correspondingly worse. Table [3] lists the derived molecular CH 3 OH and CH 3 CN column 
densities and abundances which are on average relatively low. 

The average CH 3 CN column density and abundance are 3.1 x 10 12 cm -2 and 1.2 x 10~ 10 , 
respectively. This CH 3 CN abundance compares well to the CH 3 CN abundance observed 
for example in the relatively cold and quiescent Orion extended ri dge whereas it is two 
order s of magnitude below the values obtained for the Orion hot core (Ivan Dishoeck fc Blake 
19981 ). Fu rthermore, it fits well to the chemical modeling of young high-mass star- forming 
regions by lNomura fc Millar! (12004 ): they find very low CH 3 CN column densities of the order 



10 11 — 10 12 cm -2 after only a few hundred years which then progressively increase to values 
> 10 15 cm -2 after several 10 5 years. In this picture, even the sources with detected CH3CN 
emission are as expected extremely young massive star-forming regions. 



10 



The average CH 3 OH column density and abundance are 8.2 x 10 cm and 4.3 x 10 
respectively. This abundance is orders of magnitude below the values obtained toward the 
typic al regions in Orion (e.g., 8 x 10~ 9 for the extended ridge and 2 x 10~ 7 for the hot 



core. 



van Dishoeck fc Blakdll998l ). However, it is very close to the recently rep orted CH3OH 



abundances of 6 x 10 10 in the low-mass starless cores L1498 and L1517B (ITafalla et al. 



Beuther et al. ( 20021 ) presented peak fluxes from 2D Gaussian fits to the images whereas here we use the 



exact fluxes toward the peak positions. The former are on average a bit lower. 
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20061 ). Although the statistics are still small, this similarity between CH3OH abundances 



in low-mass starless cores and high-mass IRDCs indicates that the initial CH 3 OH molecular 
abundances at early evolutionary stages do not differ significantly between low- and high- 
mass star-forming regions. 



4. Conclusions and summary 

These single-dish mm spectral multi-line observations toward a sample of 43 massive 
IRDCs revealed a range of interesting results. Since IRDC are expected to harbor sources 
in different evolutionary stages from genuine HMSCs to High-Mass cores with embedded 
low- to intermediate-mass protostars potentially destined to become massive stars, it is not 
surprising that the observational signatures do vary significantly over the sample. Because 
our spectral line selection covered tracers sensitive to different processes (outflows, dense 
gas, lower density gas, hot core chemistry) we are able to constrain various characteristics. 
In summary they are: 

• Toward 40% of the sources SiO(2-l) as a typical outflow/jet tracer was detected. 
Assuming that SiO is exclusively produced by gas phase reactions after sputtering Si- 
bearing species from dust grains, this implies that at least 40% of the IRDCs presented 
here have ongoing star formation activity (non-detection of SiO does not imply no star 
formation activity, hence the lower limit). The observed width-to- zero- velocity of the 
SiO lines varies between 2.2 and 65kms _1 . While inclination effects are unlikely to 
cause this velocity spread, we discuss other potenatial origins, e.g., embedded objects 
with different masses and/or possible evolutionary effects. 

• While we usually detect many velocity components in the low-density tracers 12 CO(2- 
1) and 13 CO(1-0), the high-density tracer H 13 CO + (1-0) is also present in all except one 
source, but in general tracing only one velocity component. Therefore, it is well suited 
to derive velocities with respect to the standard of rest and hence kinematic distances 
for IRDCs. The observed FWHM line-width of H 13 CO + (1-0) is on average a factor 
1.5 larger than that of the previously observed NH3(1,1) line. Several possibilities for 
that behavior are discussed: While the typical line-width-size relation is unlikely the 
cause for this behavior, the data are consistent with gravitational instabilities of the 
cores. However, this interpretation is not unambiguous. Since H 13 CO + has higher 
critical densities, it traces the denser central cores with likely larger internal motions. 
This motions could be caused by turbulence or by active star formation processes like 
outflows/infall, rotation. 
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• As expected CH 3 CN is only rarely detected (toward ~14% of the sources), whereas 
CH3OH is observed toward a larger fraction of the sample (~40%). Both molecules are 
only detected in their lower-energy transitions. While temperature estimates are hardly 
feasible with these data, we are able to derive the CH 3 CN and CH 3 OH abundances 
with average values of 1.2 x 10~ 10 and 4.3 x 10 -10 , respectively. The low CH3CN column 
densities are consistent with chemical models where the CH3CN abundances are very 
low at the earliest evolutionary stages and progressively increase with time. The derived 
CH3OH abundances are close to recently reported values for low- mass starless cores, 
indicating that the initial CH 3 OH abundances at the onset of star formation processes 
do not vary significantly from low- to high-mass cores. 



We like to thank Sven Thorwirth for continuous spectroscopic help. H.B. acknowledges 
financial support by the Emmy-Noether-Program of the Deutsche Forschungsgemeinschaft 
(DFG, grant BE2578). 
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Table 2:: Source parameters 



Name 






R.A. 




Dec. 




i;(H 13 CO+) 


Au(H 13 CO + ) 


Aw (SiO) b 


qc 


K 








J2000 




J2000 


kpc 


kms -1 


kms -1 


kms -1 


mjy 
f 1 1 "X 2 




IRDC18089- 


1732-3 


18 


11 45.3 


-17 


30 38 


3.6 


34.2(0.1) 


3.9(0.2) 


5.5 


147 


1.5 


IRDC18090- 


1832-2 


18 


12 02.0 


-18 


31 27 


6.6 


111.2(0.1) 


2.2(0.2) 




89 


0.9 


IRDC18102- 


1800-1 


18 


13 11.0 


-18 


00 23 


2.7 e 


21.4(0.1) 


3.4(0.1) 


42.2 


632 


6.3 


IRDC18151- 


1208-2 


18 


17 50.3 


-12 


07 54 


3.0 


29.8(0.1) 


3.0(0.1) 


65.0 


620 


6.2 


IRDC18182- 


1433-2 


18 


21 14.9 


-14 


33 06 


3.6 


40.9(0.1) 


2.2(0.1) 




144 


1.4 


IRDC18182- 


1433-3 


18 


21 17.5 


-14 


29 43 


4.6 e 


60.2(0.1) 


1.7(0.3) 


5.5 






IRDC18182- 


1433-4 


18 


21 14.0 


-14 


34 20 


3.6 


41.4(0.1) 


2.0(0.2) 








IRDC18223- 


1243-2 


18 


25 10.0 


-12 


44 00 


3.7 


45.5(0.1) 


3.4(0.1) 




156 


1.6 


IRDC18223- 


1243-3 


18 


25 08.3 


-12 


45 27 


3.7 


45.7(0.1) 


2.7(0.1) 


35.5 


289 


2.9 


IRDC18223- 


1243-4 


18 


25 06.8 


-12 


48 00 


3.7 


45.7(0.1) 


2.4(0.1) 


38.9 


77 


0.8 


IRDC18223- 


1243-5 


18 


25 12.0 


-12 


53 23 


3.7 


45.1(0.2) 


1.7(0.2) 








IRDC18223- 


1243-6* 


18 


25 16.2 


-12 


55 33 


3.7 


45.7(0.1) 


2.0(0.3) 


2.4 






IRDC18223- 


1243-6/ 


18 


25 16.2 


-12 


55 33 


3.7 


50.9(0.2) 


2.1(0.4) 








IRDC18247- 


1147-3 


18 


27 31.0 


-11 


44 46 


6.7 


59.7(0.2) 


2.2(0.6) 




110 


1.1 


IRDC18306- 


0835-3 


18 


33 32.1 


-08 


32 28 


3.8 e 


54.8(0.1) 


2.4(0.1) 




154 


1.5 


IRDC18306- 


0835-4 


18 


33 34.8 


-08 


31 20 


3.8 e 


55.2(0.1) 


2.4(0.2) 








IRDC18337- 


0743-3 


18 


36 18.2 


-07 


41 00 


4.0 


56.2(0.1) 


2.8(0.1) 


39.8 


123 


1.2 


IRDC18337- 


0743-4 


18 


36 29.9 


-07 


42 05 


4.0 


55.8(0.1) 


4.2(0.3) 








IRDC18337- 


0743-5 


18 


36 41.0 


-07 


39 56 


4.0 


58.4(0.1) 


4.1(0.2) 


16.2 






IRDC18337- 


0743-6 


18 


36 36.0 


-07 


42 17 


4.0 


56.6(0.1) 


2.5(0.2) 








IRDC18337- 


0743-7 


18 


36 19.0 


-07 


41 48 


4.0 


55.9(0.1) 


2.7(0.2) 








IRDC18348- 


0616-2 


18 


37 27.6 


-06 


14 08 


6.3 


109.7(0.1) 


2.6(0.1) 


9.9 


169 


1.7 


IRDC18348- 


0616-8 


18 


37 14.7 


-06 


17 25 


6.3 


109.3(0.2) 


4.4(0.4) 








IRDC18385- 


0512-3 


18 


41 17.4 


-05 


10 03 


3.3 e 


47.0(0.1) 


2.1(0.1) 


41.1 


77 


0.8 


IRDC18431- 


0312-3 


18 


45 45.0 


-03 


08 56 


6.6 e 


105.4(0.3) 


2.2(0.7) 








IRDC18431- 


0312-4 


18 


45 53.0 


-03 


09 01 




non det. 


non det. 








IRDC18437- 


0216-3/ 


18 


46 21.9 


-02 


12 24 


6.2 e 


97.6(0.2) 


2.1(0.4) 




96 


1.0 


IRDC18437- 


0216-3/ 


18 


46 21.9 


-02 


12 24 


7.3 e 


111.3(0.1) 


2.3(0.4) 




96 


1.0 


IRDC18437- 


0216-7/ 


18 


46 22.0 


-02 


14 10 


6.2 e 


96.1(0.1) 


2.9(0.1) 


13.0 






IRDC18437- 


0216-7/ 


18 


46 22.0 


-02 


14 10 


7.3 e 


111.2(0.7) 


1.9(1.8) 








IRDC 18440 - 


0148-2 


18 


46 31.0 


-01 


47 08 


5.6 e 


97.9(0.2) 


2.9(0.4) 








IRDC18447- 


0229-3 


18 


47 42.0 


-02 


25 12 


6.4 e 


102.1(0.2) 


2.6(0.5) 








IRDC 18447- 


0229-4 


18 


47 38.9 


-02 


28 00 


6.4 e 


99.4(0.1) 


1.9(0.3) 








IRDC18454- 


0158-1/ 


18 


48 02.1 


-01 


53 56 


3.5 e 


52.8(0.4) 


1.7(0.7) 




178 


1.8 


IRDC18454- 


0158-1/ 


18 


48 02.1 


-01 


53 56 


6.4 e 


100.2(0.1) 


2.7(0.2) 


10.9 


178 


1.8 


IRDC18454- 


0158-3/ 


18 


47 55.8 


-01 


53 34 


6.0 e 


94.3(0.1) 


3.5(0.3) 


17.3 9 


204 


2.0 


IRDC18454- 


0158-3/ 


18 


47 55.8 


-01 


53 34 


6.4 e 


98.4(0.1) 


3.3(0.3) 


17.3 s 


204 


2.0 


IRDC18454- 


0158-5 


18 


47 58.1 


-01 


56 10 


6.0 e 


94.6(0.1) 


2.7(0.1) 


17.8 


102 


1.0 


IRDC18454- 


0158-12/ 


18 


48 05.7 


-01 


53 28 


3.1 e 


46.3(0.4) 


2.5(0.8) 








IRDC18454- 


0158-12/ 


18 


48 05.7 


-01 


53 28 


6.4 e 


101.4(0.1) 


2.7(0.3) 








IRDC18460- 


0307-3 


18 


48 36.0 


-03 


03 49 


5.2 


84.6(0.1) 


1.8(0.2) 




84 


0.8 
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IRDC18460-0307-4 


18 48 46.0 


-03 04 05 


5.2 


84.4(0.1) 


2.4(0.3) 




177 


1.8 


IRDC18460-0307-5 


18 48 47.0 


-03 01 29 


5.2 


84.0(0.1) 


2.1(0.1) 




74 


0.7 


IRDC18530+0215-2 


18 55 29.0 


+02 17 43 


5.0 e 


75.9(0.1) 


2.1(0.2) 


2.2 






IRDC19175+1357-3 


19 19 52.1 


+14 01 52 


1.1 


7.7(0.2) 


2.9(0.7) 




53 


0.5 


IRDC19175+1357-4 


19 19 50.6 


+14 01 22 


1.1 


7.8(0.1) 


2.1(0.2) 




87 


0.9 


IRDC19410+2336-2 


19 43 10.2 


+23 45 04 


2.1 


21.5(0.1) 


2.3(0.1) 


36.8 h 


343 


3.4 


IRDC20081+2720-1 


20 10 13.0 


+27 28 18 


0.7 


5.7(0.1) 


1.8(0.1) 




244 


2.4 


IRDC22570+5912-3 


22 58 55.1 


+59 28 33 


5.1 


-47.5(0.1) 


2.8(0.1) 


5.5 


88 


0.9 


a Most distances are 


taken from Sridharan et al. 


(2005). 


The exceptions 


are marked. 









b Width of SiO(2-l) line down to zero intensity. 

c Here we present the 1.2 mm peak fluxes extracted directly from the peak positions. In the original paper, 

Bcuthcr et al. (2002) present peak fluxes from 2D Gaussian fits which are on average a bit lower. 
d The column density caluclations follow Bcuthcr et al. (2005a) at an average temperature of 19 K to be 

better comparable with the CH3OH and CH3CN column densities calculated at T = 18.75. Masses and 

integrated fluxes were already presented in Sridharan et al. (2005) 
e Newly derived near kinematic distances following Brand & Blitz (1993) because IRDCs at the far 

distances are unlikely. 

' Rows with the same IRDC name correspond to different H 13 CO + velocity components. 
9 This is a line-blend between both velocity components. 
" The outflow was studied in detail by Beuther et al. (2003). 
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Table 3: CH 3 CN and CH 3 OH parameters. 



Name 


Av k-0 

ft, — U 


ft, — ± 


-^CH3CN 


^CHaCN 


-^CHsOH 


^CHaOH 




kms -1 


kms -1 


10 12 cm- 2 


10 -io 


10 13 cm- 2 


10 -io 


IRDC18102-1800-1 


2.9(0.4) 


3.9(0.8) 


2.6 


0.3 


2.8 


0.4 


IRDC18151-1208-2 


3.4(0.3) 


4.5(0.5) 


5.3 


0.8 


37.2 


6.0 


IRDC18182-1433-2 










3.0 


2.1 


IRDC18223-1243-3 


2.8(0.5) 


3.1(0.5) 


2.7 


0.9 


17.6 


6.1 


IRDC18223-1243-4 










4.6 


6.0 


IRDC18337-0743-3 


3.1(0.9) 


6.2(1.4) 


4.4 


3.6 


13.1 


10.7 


IRDC18385-0512-3 










5.3 


6.9 


IRDC18437-0216-7 










2.5 




IRDC18454-0158-1 










1.7 


0.9 


IRDC18454-0158-1 










1.7 


0.9 


IRDC18454-0158-3 


76 








11.1 


5.5 


IRDC18454-0158-3 


76 


76 






11.1 


5.5 


IRDC19175+1357-4 










4.3 


5.0 


IRDC19410+2336-2 










5.0 


1.5 


IRDC20081+2720-1 


1.4(0.3) 


1.3(0.4) 


0.7 


0.3 






IRDC22570+5912-3 










2.4 


2.8 



a Line-width of the CH 3 CN k — and k — 1 components. 
b No reasonable fits because of multiple velocity components. 
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Fig. 1.— 12 CO(2-l) (grey), 13 CO(1-0) (black line) and H 13 CO+(1-0) (red line) spectra 
toward the whole sample of IRDCs. The H 13 CO + (1-0) are scaled up by the factor labeled 
in each panel for presentation purposes. 
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Fig. 2. — Figure [T] continued. 




iv(H CO (1— 0)) 



Fig. 4. — Compariso n of line- widths between H 13 CO + (1-0) (this work) and NH3(1,1) 
jSridharan et alJbpOfih . 
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Fig. 5.— SiO(2-l) detections. 
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Fig. 6. — CH 3 CN(6jt — 5^) detections. The rest frequencies of the k — 0, 1, 2 components are 
at 110383.522, 110381.404 and 110375.052 MHz, respectively. Because of the small intrinsic 
line- widths, the k — 0, 1 components are spectrally resolved which is usually not the case for 
more evolved High-Mass Protostellar Objects or Hot Molecular Cores. 
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Fig. 7. — CH 3 0(5-4) detections. The rest frequencies of the four detected lines are: 
CH 3 OH(5 1 -4 1 )E at 241879.073MHz, CH 3 OH(5 -4 )A+ at 241791.431 MHz, CH 3 OH(5_i- 
4_i)E at 241767.224 MHz, CH 3 OH(5 - 4 )E at 241700.219 MHz. The few single-channel 
spikes are artifacts caused by the correlator. 



